Electron-Positron Plasma Drop Formed by Ultra-Intense Laser Pulses 



(N 

o 

(N 



Inga Kuznetsova and Johann Rafclski 
Department of Physics, The University of Arizona, Tucson, Arizona, 85721, USA 

(Dated: September 14, 2011) 

We study the initial properties and positron annihilation within a small electron-positron plasma 
drop formed by intense laser pulse. Such QED cascade initiated plasma is, in general, far below the 
chemical (particle yield) equilibrium. We find that the available electrons and positrons equilibrate 
kinetically, yet despite relatively high particle density, the electron-positron annihilation is very 
slow, suggesting a rather long lifespan of the plasma drop. 
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I. INTRODUCTION 

Conversion of the high intensity laser pulse energy into 
a dense gas of e + , e~ electron-positron pairs is a topic of 
current theoretical and, soon, experimental interest. A 
QED cascade mechanism producing a rapid conversion 
of laser pulse energy into pairs was demonstrated in [l| 
for pulse intensity on the order of 10 24 W/cm 2 . Consid- 
ering the known reaction cross sections [2[ , subsequent to 
the electromagnetic cascade process discussed in Ref. [l| , 
photons escape the small plasma drop, while as we show 
here, the electromagnetic scattering thermalizes the mo- 
mentum distribution of this relatively dense electron- 
positron phase. We thus find a drop of '"thermal"' mo- 
mentum equilibrated, but '"chemical"' yield nonequili- 
brated electron-positron plasma with a size as small as a 
few /xm and an energy content up to a kJ. Such plasma 
will expand, and lose energy by positron annihilation. 
We obtain here the rates of energy and particle loss by 
annihilation. 

The corresponding initial local energy density is pro- 
vided by the laser field. We assume the formation of 
the plasma drop at rest in the lab frame e.g. invoking 
symmetric laser pulse collisions triggering QED cascades. 
The experimental pulse intensity parameter, defining 
plasma drop properties, is Q 



X 



a 



eE X 



in 



(1) 



where e is the electron charge, E is the laser field 
strength in the focus, A is the wavelength, and m is 
the electron (positron) mass. The discussion of physi- 
cal properties, that we present, corresponds to ciq ~ 4000 
This value will be within the range of the next genera- 
tion ultra intense pulsed lasers. For a plasma drop radius 
R = 3 /im, 2R = 3A the corresponding total plasma drop 
energy is 0(0.3) kJ. 

In the present context of plasma cooling we extended 
results of Ref. [2] to the lower density and lower tem- 
perature domain. The important theoretical refinement 
discussed here for the first time, in the context of laser 
generated low density e~e + plasma, is the consideration 
of the plasmon screening depending on plasma tempera- 
ture and density. We also extend our earlier considera- 



tions to the nonrelativistic regime T < m as required in 
the study of the plasma expansion and freeze-out process. 

Under the experimental conditions we consider here, 
all photons produced will escape from the small drop of 
low density plasma of electrons and positrons without 
much, if any, scattering. However, even far from the 
chemical equilibrium density of the particle pair yield, 
it is possible for the produced electrons and positrons 
to equilibrate thermally by means of M0ller and Bhabha 
scattering, 



e ± + e ± O e ± + e ± , 
e ± + e T <^ e ± + e T . 



(2) 
(3) 



forming an electron-positron plasma drop: wWhen the 
drop size R exceeds the scattering length L cc , 



R> L c 



(4) 



multiple scattering processes can occur, allowing ki- 
netic '"thermal"' equilibration. We therefore study 
positron annihilation loss processes assuming the Fermi- 
Boltzmann energy distribution of available particles. We 
solve kinetic population equations and evaluate the frac- 
tion of particles in plasma which can annihilate during 
the plasma life span. 

There are two paths to positron annihilation, the direct 
in-flight pair annihilation, 



e ± + 



7 + 7 



(5) 



and in-flight bound state positronium p s formation, 



e ± + e T -> 7 + Pn Pn ^717, n = 2, 3 (6) 

which is followed ultimately by annihilation. The anni- 
hilation life span of positronium for spin is rp 2 = 0.12 
ns, while for spin 1 it is Tp3 = 140 ns. However, the 
positronium formation cross section only competes with 
the in-flight annihilation cross section for temperatures 
below T w 60 eV Qj], and at that point, the expansion 
dilution will, in general, slow these processes down con- 
siderably. 

In Sec. [XT] we present cross sections for M0ller and 
Bhabha scattering including in the plasmon screening ef- 
fects. We compare the resulting pair annihilation cross 
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section with positronium formation. In Sec. IIIII we 
present numerical results for the M0ller and Bhabha scat- 
tering mean free path and also annihilation relaxation 
time. We discuss conditions for plasma drop to be ther- 
mally equilibrated. In Sec. IIVI we evaluate our results 
and present conclusions. 



where V = AirR 3 /3 is the volume and g e ^ = 2 is the spin 
degeneracy. When the e,e-pair yield is far below chem- 
ical equilibrium, that is, T -C 1, the effects of quantum 
statistics are, in general, less significant and the Boltz- 
mann limit is often equally precise, 

/ e , e - -> T e - e e~ E ' T . (11) 



II. e+, e - PLASMA REACTION RATES 
A. Scattering rates 

1. Particle Density 

We consider the case of a small nonopaque expand- 
ing electron-positron plasma drop. The drop stays ther- 
mally equilibrated by scattering processes. The electron 
(positron) multiplicity Ni (i — e + ,e~) is thus in ther- 
mal (momentum distribution) but not in chemical (yield 
distribution) equilibrium. 

It has been shown Q that in order to maximize the 
entropy at fixed particle number the appropriate maxi- 
mum entropy distribution is the usual Fermi-Dirac / ej g 
distribution accompanied by a phase space occupancy 
parameter T, 



fe 



1 



(7) 



T(t) describes the pair density and is, in general, a func- 
tion of time, and it is the same for both particles and 
antiparticles. This is in contradistinction to the chem- 
ical potential v which changes sign, v s = —v e compar- 
ing particles and antiparticles. The chemical potential v 
regulates the abundance difference between particles and 
antiparticles and thus, in general, is only weakly depen- 
dent on time. A system with T = 1 for all particles is 
in chemical equilibrium, and we refer to particle density 
with T = 1 as a chemical equilibrium density. 

Note that the Lorentz-invariant exponents involve the 
scalar product of the particle four-momentum with 
the local four- vector of velocity it' 1 , where w M describes 
the local collective flow of matter, as expected for an 
unconfined plasma drop. The thermal properties is, T, T 
are defined in the local rest frame. In the absence of local 
matter flow the local rest frame is the laboratory frame, 



2. Plasmon mass and screening length 

To avoid Coulomb singularity in reaction matrix ele- 
ments we introduce the plasmon mass, induced by the 
plasma screening effect, following the example of gluon 
dynamics in quark-gluon plasma @. The plasmon mass 
is0 



je+ + fe- ( 1 _J^_\ 



3E%J (2tt) 



(12) 



a = e 2 /47r = 1/137.036 is the fine-structure constant. 
For nonrelativistic temperatures T « m e , m 1 goes to 
the classical plasma frequency, and a simple limit also 
emerges for relativistic temperatures with T = 1, 



47ry / 2an e /m e T < m e 



V4toT/3 



T > m e , T = 1 



(13) 



The corresponding screening length, the Debye radius, is 

vt 



and the mean thermal particle velocity vt is 



Vt 



(14) 



(15) 



since f e + — f e - — f. We show in the Fig. Q] the 
electron (positron) screening length and the mass of the 
plasmon as a function of T. The plasmon mass is in- 
creasing towards the small temperatures and is asymp- 
totically constant, similar to the behavior of the plasma 
density. The screening length is otherwise decreasing to- 
wards the small temperatures (inverse proportional to 



and vt oc v T) in our range of temperature. 



We thus have 



1 



/e ' 5 T L r!e B / T + 1 

The yields of particles are 



N e ,e = n e ,eV = 9e,eV 



d 3 p 
(2^)3 



(8) 



(9) 



(10) 



3. Boltzmann limit 

We are interested in experimental conditions under 
which the number of pairs produced is large compared to 
the residual electron density originating in matter. Fur- 
thermore we will deal with conditions (T c < 1 or/and 
T < m MeV) which allow us to use the Boltzmann ap- 
proximation. Then, we have 



Tie 



sinh^/T) < 1. 



(16) 
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FIG. 1: Upper panel: electron (positron) screening length as 
a function of plasma temperature. Lower panel: mass of the 
plasmon as a function of T. 



In what follows we will set v — 0, and consider elsewhere 
the case for very low density degenerate plasma, where 
the chemical potential may become important. We thus 
have T ejg = T. In the relativistic Boltzmann (classical) 
limit the plasma density and energy density are 



T e g e T 3 



T 



2tt 2 
3, 9e T 4 
2 2ir 2 



x 2 K 2 (x), 



(17) 



x 2 K 2 (x) + -x 3 K 1 {x) ) , (18) 



where Ki(x) is a Bessel function, x = m/T. 

4- Electron (positron) scattering rates 

In the evaluation of the matrix element we use Man- 
delstam variables: s, u, and t. In the case of M0ller 
scattering 

s = (Pi+P2) 2 ; u= (p 3 -p 2 ) 2 ; t = (p 3 - px) 2 ; (19) 

and s + u + t = m\ + m\ + m\ + m\ . 

The M0ller scattering matrix element is [1, H, @ , 



\Mr± P ±\ 2 = 2W 



,s 2 + u 2 + 8m 2 (i - m 2 ) 
2(i-m 2 ) 2 



+ t 2 + 8m 2 (u-m 2 ) (s ~ 2m 2 ) (s - 6m 2 ) 



2{u — m 2 ) 2 



(t — m 2 )(u-m 2 ) [ 



In the case of Bhabha scattering we have 

s = (P3-P2) 2 ; U = { Pl +p 2 ) 2 - t = {p 3 -pi) 2 , (21) 

see diagrams in 0] . The matrix element does not change 
in terms of variables Pi,p 2 ,P3, when it is written in terms 
of variables s, u, t we need to cross u and s in the Ivfoller 
scattering matrix element [see Eq. (|20|) ]. 



|M e±eT (M, U )| 2 = |M e±e± (u,M)| 2 ; (22) 



thus we find 



|M e±eT | 2 = 2 b 7r 2 a 



(! ., , s 2 + u 2 + 8m 2 (t - m 2 ) 



2(t~m 2 ) 2 



+ t 2 + 8m 2 (s - m 2 ) (u - 2m 2 ) (u - 6m 2 ) 

(t - m 2 )(s - m 2 ) 



2{s-m 2 ) 2 



>(23) 



For M0ller and Bhabha scattering the cross section 
<7 e e(s) can be obtained by averaging the matrix element 
over the t variable: 



<Jee{s) 



\6-k(s — Am 2 ) 2 



dt\M ee \ 2 , (24) 



where t min = — (s — 4m 2 ), t max = in both cases @. 
Similar evaluations were done for heavy quarks produc- 
tion pjj. 
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For M0ller and Bhabha cross sections we obtain in plasma, keeping m 7 , 



CT e ± e ±^ e ± e ±(s) 



1 



dt\M t 



Ana 2 



167r(s-4m 2 ) 2 7_ (s _ 4m2) ee (s - 4m 2 ) 

87ra 2 ( (s — 2m 2 ) (s — 6m 2 ) 



s 2 + 8m 2 (m 2 — to 2 ) + (s + m 2 — 4m 2 ) 
(s + to 2 — 4m 2 ) m 2 



i s ~\~ nr 2 I In — 

(s — Am 2 ) 2 \_ [s — Am 2 + 2m 2 ) ' 1 J s — Am 2 + to 2 



(25) 



Cr e ±e=F« e ±eT(s) 



i r 

167r(s-4m 2 ) 2 ]_, 



— (s— 4m 2 ) 



tft|Afe± oT | 



27TQ! 2 



(s — 4m 2 ) 



s 2 + 8m 2 (m 2 - m 2 ) + (s + m 2 - 4m 2 ) 2 ( ^ 8 ((s - 4m 2 ) 2 + m 2 (s - m 2 )) + 3s + 2m 2 + 4m 2 



4m 2 ) 



3(s - m 2 ) 2 



(s - m 2 ) 



+2 



(m 2 + s) 2 — 4m 4 + (s 2 — to 4 ) 



7 ' In 



s — Am 2 + to 2 



r 



(26) 



B. e + e Annihilation 



Master equation and annihilation time constant 



The master population equation reads 

1 ^ = — TpTgWai 



1/ eft ' E ' E " ra1, (27) 

We have made explicit the dependence of evolution of the 
particle (pair) multiplicity in thin plasma on the prevail- 
ing density showing the factor T e T g . 

A simplified form of the master equation (up to dilu- 
tion by volume expansion, to be considered elsewhere) is 
easily obtained, 

1 dT e 1 T g 

T r U e T in ' 

1 e U* ' ann 1 e 

introducing the annihilation relaxation time rf nn 
_ rfn e /dT c 

'ann TinW ' y ' 

L e yv ann 

and similarly for rf nn . In our case T e ~ T g and we see 
that 

TiP /"* <ft' 

^r= / 4h*0 (so) 

1 e JO r ann 

We can write a similar master equation for the plasma 
drop energy loss, 

1 (t ptot 

^^ = -T e T e -^ n , (31) 

where E tot is the total energy of the plasma drop. The 
relaxation time of energy loss is 

E de/dT e 



T in W E ' 



(32) 



where e is the plasma energy density and T™ is the initial 
electron (positron) phase space occupancy. 



2. Annihilation rate in flight 

When electrons collide with positrons, they can anni- 
hilate. We consider here the dominant in flight annihi- 
lation process into two photons. The invariant rate of 
annihilation per unit of volume and time e + e — > 7 + 7 
is (3 + 4 -> 1 + 2) 



97 



2(27r) f 



d 3 P j r d 3 P 7 2 r d 3 pi r d 3 P : 



2Ej 



2El 



IE* 



2£f 



^ 4 (P^ +Fl - Vl - Pi) E I (PlP2 iMj^eelptpl) 



spin 



xe -(rf+^)/T /e(? ,e )T -l /g(p e )T -l_ 



(33) 



Here (pip? |Af 77 ^. e g| P3P4) is the annihilation quantum 
matrix element which we will consider to lowest order 
in a, g e is electron-positron degeneracy, and factor 1/2 
is due to the indistinguishability of the final state pho- 
tons. We used this method to describe the electron- 
positron pair annihilation in [2j], adapting it from work 
on strangeness production in quark-gluon plasma [TTI — 
[H]. In the last line of Eq.flHHJl we introduce T~ 1 TZ; 1 to 
compensate for the factor T e Upsilon g seen in Eq. (|27l) . 

The invariant rate Eq. ([33)) relates to the electron- 
positron pair annihilation cross section [l5[ ; in the Boltz- 
mann limit we have 



Want, = 



32tt 4 



ds\fs(s — Am 2 ) a 



ee ^ 11 (s)K 1 (^/s/T). 

(34) 
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Here the annihilation cross section is Q 



^"ee- ^77 (^) 



27ra 2 (s 2 + 4to 2 s - 8m 4 ) 
s 2 (s — Am 2 ) 



y/s + Vs- Am 2 (s + Am 2 )\fs 2 - 4m 2 s \ 
n ^T7^1^f ~ ( S 2 + 4m 2 s _ 8m 4) I -i 35 ) 



5. Energy loss 

Once in-flight e + e annihilation occurs, the produced 
photons escape the small plasma volume. An analogous 
expression to Eq. ([55)1 describes the energy loss rate due 
to pair annihilation, 



W E = — ^ 

' ' n m <n 



g 2 f^pj fd 3 pl f d 3 pl f d 3 p] 



2(2tt) 8 J 2Ej J 2E] J 2E% J 2E% 

xs* (jpI+pI-pI-pDYjMpI \M ee ^\ptpi)t 



x (El + Et)f e (pt)f e (pt)T- 2 e^+P^ T , (36) 

We now obtain a relation analogous to Eq. (|34|l . Con- 
sider the integral [l5[ leading to Eq. (f3~4"|) , 



d 4 pe-^ u S ( P 2 -s) = jVtKtfVs), (37) 



where u = (1,0) in the laboratory frame. Instead, we 
now need to use 



J d 4 pp-ue-^ u S Q (p 2 -s) = -^JlyTsK^yTs). 



We use d[K\{x)/x]/dx = —K 2 (x)/x to obtain 



(38) 



W. 



g 2 T 



dss(s — Am 2 ) a ee ^ 11 {s)K2{yfs/T). 



~ 32^ 4 ' ~"" y " '""-r^ 



4- Positronium formation 



(39) 



The cross section for radiative positronium (ee) forma- 
tion, e~ + e + O 7 + (ee) is 



2 12 7r 2 w / e \ e 



3pm 2 V 1 + £ 2 / 1 - e ~ 27Ti 



-,— 4^arccot^ 



5p 2 



(40) 

where £ = am/2p and the photon energy u) is defined by 
the conservation law 



is valid while £ < 1. This condition is satisfied up to 
temperatures on the order of 10 eV. 

We did not consider in detail the influence of plasma 
screening on positronium formation, a topic which invites 
further work in view of currently available results. It was 
found in [llj that the plasma screening and collective 
effects significantly reduce the radiative recombination 
cross section in non ideal plasma. The screening effect 
for positronium formation should be similar to result for 
free electron radiative recombination with ions in non- 
ideal classical plasmas. However, in positron - hydrogen 
plasma the Debye screening can result in a large increase 
of the positronium formation cross section at incident 
positron energy 20-100 eV |18|. 



III. RESULTS FOR LASER FORMED PLASMA 

A. Parameters for thermal plasma drop 

We assume here that the total energy E of (colliding) 
laser pulses converts in the initial volume V to the e + e~- 
plasma drop energy. The initial energy density e = E/V 
is obtained from Eq.([T|) and is characterized by ao and A, 



1 



2 _ _t_ / a m y 



Att An V eA I 
The phase space occupancy of the plasma drop is 



T = 1 (aomy 
e AireoiT) \ eA J 



(42) 



(43) 



where we introduced the chemical equilibrium energy 
density eo = e|x e =i) Eq. (fT5|l . Then the total energy of 
plasma, E, is defined by the plasma drop radius R for a 
given parameter ao and wavelength A. The initial plasma 
size is expected to be close to the wavelength. We take 
the wavelength R = 3 A/2 for all cases considered below. 

In Fig. [5] we show the phase space occupancy T e from 
Eq. (|4"31 (upper panel) and the corresponding plasma den- 
sity n p i (lower panel). The solid (blue) line shows the 
actual chemical nonequilibrium values. For comparison 
the chemical equilibrium results are shown by the dashed 
(green) line. We note that for T » 0.06 MeV the fully 
equilibrated yield is much greater than what we can make 
using a near future high intensity laser. However, the 
density of particles in plasma which we achieve is very 
high. 

At T << m, when plasma becomes nonrelativistic the 
energy /particle — > mc 2 is a constant and does not depend 
much on the plasma temperature. Hence, the plasma 
particle density goes for T — >• to a constant for a given 
energy and plasma drop size, 



u> H = p 2 /m + a 2 m/A. 

Am 



(41) 



p is the electron (positron) momentum in the center of 
mass reference frame, p = Vs — Am 2 /2. Equation (|40[) 



n p i 



(44) 



and temperature cannot be determined considering a 
given available energy constraint. 
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FIG. 2: Upper panel: electron (positron) phase space occu- 
pancy T e as a function of T for ao = 4000 and R — 3 /im 
(solid blue line). Lower panel: plasma density corresponding 
to the phase space occupancy on the upper panel (solid blue 
line) and equilibrium density T = 1 (dashed green line) as a 
function of T. 



In a system where particle (pairs) can be produced but 
energy is fixed the entropy density reaches maximum at 
T = 1. We show the entropy density of electron - positron 
plasma, 

S = J S ((/e " 1} ln(1 " /e) " /e ln(/e)) ' (45) 

at E = 0.3 kJ and R = 3 fim as a function of temperature 
in Fig. [3] As expected, the maximum of the entropy den- 
sity is at the temperature, T = 0.06 MeV, where phase 
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FIG. 3: The entropy density of electron-positron plasma with 
ao = 4000 and R = 3 am as a function of temperature. 



space occupancy of electron and positron T e = 1. How- 
ever, the maximum is very flat. Note that there is much 
less entropy density when the system is formed at rela- 
tively high temperature. This is because there are fewer 
particle pairs and, for a relativistic gas, the entropy per 
particle is near S/N ~ 4. For far off equilibrium low den- 
sity systems the expansion of the volume is thus accom- 
panied by reactions that tend to chemically equilibrate 
the system and move it towards chemical equilibrium. 



B. Electron and positron scattering 

The formation of electron-positron plasma is further 
subject to the opacity condition Eq.Q. To check if 
this condition is satisfied we extend our earlier consid- 
erations [||, now introducing plasmon mass, Eq. (fT!?|) . in 
a domain of mild relativistic and nonrelativistic temper- 
atures. 

The electron (positron) mean free path follows from 



n e 



(46) 



where for the scattering rate W cc we use an equation 
similar to Eq. (|3~4"|) (since the final state does not have two 
identical bosons, the normalization factor is different): 



Wee = 



and 



9 2 T 
32tt 4 



4m a 



dsVs(s-4m 2 )a ee (s)ifi(Vs/T), (47) 



(48) 
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Bhabha scatterings, Eq. (|25j) and Eq. (|26|) . are 
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FIG. 4: Electron (positron) scattering length at a given 
plasma radius and energy as a function of T. 



In Fig. S] we show the electron (positron) scattering 
length L CCl Eq. (|46[) . at a given plasma radius R = 3 (im 
and energy 0.3 kJ (a = 4000) as a function of plasma 
temperature T. T varies for every value of T, as we see 
in Fig. [2] Since T << 1 the scattering length can be 
evaluated in the Boltzmann limit in practically the entire 
temperature range of interest, including T > m. We 
also show (dashed green line), for comparison, the case 
T c = 1, which means that we allow the density to go up 
significantly and the small difference we see in figure 2] 
for high T is due to quantum gas properties. 

At relativistic temperatures T ~ 1 MeV our present re- 
sult is in agreement with scattering rates evaluated with 
plasmon mass taken in the limit of ultrarelativistic tem- 
peratures in [2| with an accuracy of few percent. 

For the constant plasma drop energy scattering length 
Lee bas a maximum at T w m . In the whole temperature 
range the plasmon mass is small and the first term in 
Eq. (|25p and Eq. (|2"B]) is dominant, resulting in the cross 
section for electron or positron scattering, 



(Tee ocm 7 an 



(49) 



In the range where condition (l4T)|) is valid, the electron 
(positron) mean free path does not depend on density or 
T e . When the mean free path is increasing with decreas- 
ing density, this is compensated by a larger cross section 
because of a smaller plasma screening effect or smaller 
to 7 . For the entire T range, the scattering length scale is 
a tiny fraction of the plasma size. 

In the temperature range T < to the contribution of 
Ap 2 = s — 4m 2 is much smaller than m 2 and much larger 
than to 2 ; the approximate cross sections for M0ller and 



(T e ± e ±^ e ± e ±(s) = 2er e ± eT ^ e±e ±(s) 



6Ana to 



(s — Am 2 ) 2 to. 2 

(50) 

One can also consider a Rutherford-type differential 
cross section for M0ller scattering Q 



dcos( 



Ap 4 



-cosec 8/2. 



(51) 



We checked by integrating Eq. (|5"Tj) numerically that 
Eq. (|50| corresponds to the total cross section from the 
integrated Eq. (|51|) with a cutoff angle 9 m in — rrij/m. 

We found from the results presented in Fig. 2] that 
condition [4] is satisfied for whole temperature range con- 
sidered. We conclude that the electron-positron plasma 
drop can stay thermally equilibrated at relatively low 
densities when T << 1 and/or the temperature T « to: 
the electron-positron mean free path decreases when the 
temperature decreases below the electron mass because 
of the factor s — Am 2 — Ap 2 in the denominator of cross 
section Eq. (|50| . At a temperature higher than m the 
other terms begin to contribute to the cross sections, 
Eq. (|2"5)) and Eq. (f2"o| . The electron-positron mean free 
path decreases again. 

The cross section, Eq. (|50[) , is valid in the temperature 
range 



T cr = — <T < m, 

TO 



(52) 



The reader should keep in mind that the present consid- 
erations do not automatically apply to the case of a de- 
generate electron-positron gas (high density or/and low 
temperature), where we should extend the investigation 
of collective plasmon dynamics in order to obtain a valid 
estimate of the electron-positron scattering cross section. 



C. Annihilation 

1. Plasmons 

While screening and plasma oscillations impact the 
scattering processes, this is not the case for our domain 
in regard to the annihilation process. There are several 
processes to consider: 

1. The electron (positron) thermal mass correction, 
which is on the order of magnitude of m 7 . However, 
and this correction is small. 

Plasmon f-> e + e~, if the reaction threshold is ex- 
ceeded, to 7 > 2to 6 Q. This can only happen at 
ultrarelativistic temperatures. In the case consid- 
ered here with constant plasma energy, m 2 oc T^ 1 
(see Fig. [IJ, the threshold condition cannot be sat- 
isfied. 
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3. The hard photons from annihilation (fc w m) are 
rescattering on plasmons. The condition where 
screening has noticeable effect on the photon prop- 
agation is [m 

kr D < 1, (53) 

where k is the photon wave number and m is the 
Debye radius Eq.flTJJ. This condition is equivalent 
to the condition T < T cr , Eq. (|52|) . We do not con- 
sider here such low temperature plasma. 



2. Annihilation life span 

We determine, using the perturbative QED reaction 
rate, the annihilation rate of plasma under the condi- 
tions considered in the previous subsections. We assume 
that the plasma drop formation life span is on the or- 
der of magnitude of the laser pulse duration, 10 fs, and 
this is the stage at which the density of pairs and thus 
annihilation should have the largest rate; however, this 
is not the case, since, as T increases, the pair density 
drops, given the constant initial total energy, and thus 
the annihilation relaxation time increases. 

In Fig. [5] we show relaxation times r for particle num- 
ber annihilation T ann (thick lines) and energy loss r^ nn 
(thin lines) for plasma at ao — 4000, E = 0.3 kJ (solid 
blue lines) and ao = 8000, E = 1.2 kJ (dashed green 
lines) as a function of temperature. The values of r are 
indeed largest for initial highest temperatures and there 
is a shallow minimum at T « 0.065 MeV. At T < 0.065 
MeV the pair density is approximately constant but par- 
ticle temperature decrease results to increase of annihi- 
lation relaxation time. The fastest annihilation occurs 
here because we have at this low temperature the high- 
est mobility of particles at high density. 

We recognize that the fraction of annihilations is very 
small initially, we obtain from Eq. (|27|) 

N &nn /N » T 2 W ann — « — • (54) 

Another way to look at the conditions fo annihilation is 
to note that the relaxation time is inversely proportional 
to T c . Then from Eq.([4"2)) we have 

A 2 

T ann ^ ~ o> (55) 
a 

which explains the dependence on ao (see Fig. [5]). 

We see in Fig. [5] that the energy loss relaxation time 
T ann becomes very close to r ann for T < m, since the 
energy of the plasma drop changes mostly because of the 
pair mass disappearance and the resulting decrease in 
plasma mass . At T > 2m, the energy loss relaxation 
time is, as expected, above the annihilation relaxation 
time. This happens since there is a preference for slower 




FIG. 5: Time constant for particle annihilation (thick lines) 
and energy loss (thin lines) at ao = 4000, E = 0.3 kJ (solid 
blue lines) and ao = 8000, E = 1.2 kJ (dashed green lines) as 
a function of plasma drop temperature. 

particles to annihilate, and thus on average, in the ther- 
mal bath few particles of higher energy lost and annihi- 
lation leads to a slight increase of the ambient plasma 
temperature. 

Our result seen in Fig. [5] implies that the annihilation 
process, even at the highest initial density, is relatively 
slow compared to other dynamical effects controlling the 
plasma drop: the plasma drop must live t » T ann to 
have most positrons in the plasma annihilated. This time 
is much longer than the pulse duration, 10 fs; indeed, 
it is on the scale of nano seconds. There is, further- 
more, the kinetic expansion leading to further dilution 
of the plasma - when the plasma drop expands with 
time, the density decreases and the annihilation relax- 
ation becomes even longer. Most, if practically not all 
of the 3 x 10~ 4 — 10 -5 annihilation events originate in 
the densest plasma stage during laser pulse, and a reliable 
prediction of the total annihilation yield requires detailed 
control of the kinetic processes in the initial state of the 
plasma as well as a precise understanding of the plasma 
drop expansion dynamics, which further reduces the an- 
nihilation rate, ultimately leading to a cloud of streaming 
electrons and positrons. 

3. In-flight annihilation compared to positronium formation 

In Fig.[6]we compare the nonrelativistic limit of the an- 
nihilation in-flight cross section (dashed line) to the cross 
section for radiative positronium (ee) formation (solid 
blue line) as a function of electron(positron) kinetic en- 
ergy in the center of mass frame -Ekin = (s — Am 2 ) /8m. 
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FIG. 6: Radiative positronium formation (solid line) and di- 
rect annihilation cross sections (dashed line) as a functions of 
electron (positron) kinetic energy in the two particle center 
of momentum frame. 

The cross sections intersect at Ek in w 150 eV. This cor- 
responds to the crossover temperature obtained in 
T e ~ 60 eV. Thus the direct annihilation dominates down 
to this low temperature, and our prior results apply for 
T > T e . For T < T e we have significant positronium for- 
mation only if we reach this condition without much of 
expansion, which is not part of our present study. 

IV. CONCLUSIONS 

The key result of this study is that high intensity QED 
cascading leads to an electron-positron drop which does 
not annihilate but thermally equilibrates. In this plasma 
drop electron - positron pairs are thermalizcd by M0ller 
[Eq.©] and Bhabha [Eq.Q] scattering, and they anni- 
hilate very slowly; see Fig. 

We found that in the Boltzmann limit the electron 
and positron scattering length nearly does not depend on 
plasma density in considered temperature range due to 



collective plasmon effects. The cross section decrease at 
lower density is compensated by plasmon charge screen- 
ing in the less dense plasma. As a result electron-positron 
plasma can be thermally equilibrated at the density and 
temperature range considered, far below the chemical 
equilibrium of the pair yield, T = 1. The plasma drop 
size allows very many scattering processes; we did not 
find any restriction on the minimum plasma drop energy 
and/or maximum drop size by considering the opaque- 
ness condition Eq.Q for electron (positron) scattering. 

We calculated, as an example, the annihilation relax- 
ation time for an internal plasma drop energy of 0.3 - 
1.2 kJ and radius 3 fjxa. Because of the relatively low 
density the annihilation relaxation time is much longer 
than the pulse duration, which is ~ 10 fs. We found that 
in-flight annihilation is fastest at T = 0.065 MeV, yet 
still relatively slow. The radiative positronium produc- 
tion process exceeds the in-flight annihilation at a much 
lower temperature, 60 eV, leading perhaps to the for- 
mation of positronium in the late stages of the drop. If 
such a low temperature is reached without drastic expan- 
sion dilution, very many positroniums can be formed, and 
positronium formation prolongs the life span of positrons, 
though the nature of the plasma drop is now different. 

The experimental conditions will determine at what 
temperature and, more importantly for the following ar- 
gument, rapidity, relative to the laboratory frame of ref- 
erence, the electron-positron drop will be formed (20l[2l|. 
Multipulse arrangements can be easily obtained, result- 
ing in the plasma drop being formed at high rapidity. 
The greater the rapidity, the greater the effect of the 
time dilation that prolongs the life span of the plasma 
drop, as seen in the laboratory. We recognized, in this 
work, the relative stability against annihilation evaluated 
in the intrinsic rest frame of the drop. Therefore, it ap- 
pears possible to create, using high density lasers, a qua- 
sistable matter-antimatter plasma drop capable of trav- 
eling macroscopic distances before dissipating into a low 
density cloud of particles. 
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